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A self-consistent theory of shot noise in ballistic two-terminal conductors under the action of long- 
range Coulomb correlations is presented. Analytical formulas for the electron distribution function 
and its fluctuation along the conductor, which account for the Coulomb correlations, have been 
derived. Based upon these formulas, the current-noise reduction factor has been obtained for biases 
ranging from thermal to shot-noise limits as dependent on two parameters: the ratio between the 
length of the sample and the Debye screening length A = d/Lo and the applied voltage qU/ksT. 
The difference with the formulas for a vacuum diode is discussed. 

PACS numbers: 73.50.Td, 72.70.-fm 

I. INTRODUCTION 

Recently, significant attention has been focused on 
the study of nonequilibrium fluctuations of current (shot 
noise) in mesoscopic conductors.EJ The term "shot noise," 
appearing originally in the context of pure.ballistic elec- 
tron transmission in vacuum-tube devices ^ has acquired 
nowadays a much broader usage and refers to different 
mesoscopic structures, including diffusive conductors and 
resonant-tunneling devices, where the carrier flow ex- 
hibits njonequilibrium noise proportional to the electric 
current. El 

A matter of particular interest is the significance of 
long-raiige Coulomb correlations in the noise-reduction 
effect .Bu Coulomb interactions may keep nearby electrons 
apart and more regularly spaced rather than strictly at 
random, which leads to the noise reduction, as pointed 
out by Landauer.B This effect occurs in different physi- 
cal situations. Among them are charge- limited ballistic 
transport, resonant tunneling, single-electron tunneling, 
etc. For the ballistic conductors an electrostatic potential 
barrier is formed near an injecting contact. The barrier 
fluctuates synchronously with random electron passages 
through it, which leads to noise reduction, as evidenced 
recently by Monte Carlo simulations for semiconductor 
ballistic diodesJj In this way, an incoming Poissonian 
flow is converted into an outgoing sub-Poissoniap, flow, 
exhibiting a motional electron-number squeezingB This 
effect is similar to iihat leading to shot-noise suppression 
in vacuum diodes.Ll^EI Under the resonant tunneling ef- 
fect, a built-in charge inside a quantum well affects the 
position of the resonant level and prevents the incoming 
carriers from passing through the well, thereby resiil 



ing in carrier correlation and-shot-noise reductionll2rt3 
in a certain range of biases.EjLj The Coulomb correla- 
tions in these systems act under the coherent as well 
as under the sequential tunneling regime of the carrier 
transport. The carrier correlations reach their extreme 
form of the Coulomb blockade of the electron transfer 



under the single-electron tunneling effect- leading to the 
noise reduction-studied theoreticallyli3"c3 and observed 
in experiment .EiJ 

All the above-mentioned cases have the common fea- 
tures which are necessary for the Coulomb regulation ef- 
fect and shot-noise reduction in the whole frequency spec- 
trum to occur: (i) the existence of a potential barrier in- 
side a device or at the interface with an injecting electron 
reservoir, which controls the current; (ii) the dependence 
of the barrier height and/or carrier transmission on the 
current. If no barrier is present, no shot-noise reduction 
at low frequencies due to Coulomb repulsion is expected. 
At high frequencies, however, the noise level may also be 
affected by Coulomb coctelations due to screening in an 
external environment E3'E^ 

The potential barrier, which controls the current, 
appears in an ordinary situation of the space-charge- 
limited transport. For ballistic nondegenerate conduc- 
tors this cajSehas been treated recently by Monte Carlo 
simulationsOij and attracted some attention in Ref . ^ for 
a degenerate case. For the case of diffusive nondegener- 
ate conductors, studied by the Monte Carlo technique in 
Ref. the self-consistent kinetic theory of noise, which 
takes into account Coulomb correlations, has been devel- 
oped recently in Refs. and ^ A similar kinetic theory 
for the ballistic case is lacking. 

It is the aim of this paper to address the problem of 
Coulomb correlations in ballistic conductors and present 
for the first time a self-consistent theory of shot noise in 
these conductors by solving analytically the kinetic equa- 
tion coupled self-consistently with the Poisson equation. 
It is important to compare the present noise theory for 
a semiconductor ballistic dwde with that for a vacuum 
diode developed long ago.Eru The main advance for the 
latter has been done in the celebrated paper by North 
published in 1940, where he derived an asymptotic for- 
mula for the pCurrent-noise spectral density at the high 
voltage limit Monte Carlo aiinulations of noise in vac- 
uum diodes are also availableES^ It should be stressed. 
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however, that despite of the similarity of the underly- 
ing physics (in both cases the nondegenerate Boltzmann 
electron gas without collisions in the electrostatic field 
is under consideration), the case of the semiconductor 
diode differs by several features: (i) due to a two-terminal 
geometry, there are two opposing currents instead of a 
single current, which results in different current-voltage 
characteristics at low and moderate biases;El (ii) the bal- 
listic transport regime is limited by the presence of dis- 
order, impurities, etc. Even in a pure and perfect solid, 
carriers may interact with a lattice (phonons), which at 
high biases becomes significant and breaks down the bal- 
listic regime. This makes it practically impossible to at- 
tain in solids the regime where the known formulas for 
vacuum electronics, such as the Child law for I-V char- 
acteristics or North's asymptotic formula for the noise, 
may be applied. This issue will be addressed in the pa- 
per, using the derived formulas and considering them in a 
full range of biases. Finally, we suggest an electron spec- 
troscopy experiment to make the Coulomb correlations 
effect observable. The possibility of such an experiment 
is based on recent advances in nanoscale iaibrication tech- 
niques and shot noise measurements.E3~L3 

The paper is organized as follows. In Sec. II we de- 
scribe the semiconductor ballistic structure and discuss 
the main assumptions concerning underlying physics. In 
particular, the validity of the one-dimensional plane ge- 
ometry approximation for the fluctuation problem is ad- 
dressed. In Sec. Ill we introduce the basic equations that 
describe the space-charge-limited semiclassical transport: 
the coUisionless kinetic equation coupled self-consistently 
with the Poisson equation. The steady-state problem is 
solved in Sec. IV, and the results are compared with the 
Monte Carlo simulations. In Sec. V we solve analytically 
the fluctuation problem and derive the formula for the 
current-noise spectral density that covers the range of 
biases from thermal to the shot-noise limits. The results 
for the noise-reduction factor are compared with Monte 
Carlo simulations and North's asymptotic formula for 
vacuum diodes. The contributions of different electron 
energy groups to the noise are found, and the correlations 
in energies for the electrons collected at the receiving 
contact are discussed. Finally, Sec. VI summarizes the 
main contributions of the paper, and in the Appendix we 
present mathematical details concerning the derivation of 
the fluctuations of the electron distribution function in 
the self-consistent electric field. 



II. THE PHYSICAL MODEL 



Before proceeding with a discussion of the problem, 
we will specify the structure under consideration and 
the main assumptions concerning the underlying physics. 
Consider a two-terminal semiconductor ballistic sample 
with plane parallel contacts a.t X — and X — d (see 
Fig. |l). The contacts, which we denote by L and R (left 



and right), are assumed to be heavily doped semiconduc- 
tors with a contact electron density much higher than 
that in the sample. 




FIG. 1. Schematic band-energy diagram for a n-i-n bal- 
listic diode under a space-charge-limited conduction. Two 
different types of the contacts are shown: (a) homojunctions; 
(b) heterojunctions. Shadowed regions illustrate the energy 
distribution function of electrons at the contact-sample inter- 
faces. 



The structure may then be considered as a n-i-n diode 
operating under a space-charge-limited current regime in 
which the current is determined by a charge injection 
from the conlyayts rather than by intrinsic carriers of the 
active regionEII Two different types of the contacts may 
be considered depending on whether the contact and the 
sample are fabricated of the same or different material. 
For the former case the diode is composed of two homo- 
junctions [Fig. 0(a)], while for the latter, it is composed 
of two heterojunctions with a jump of the conductance 
band Sc at the contact-sample interface [Fig. |l|(b)]. The 
underlying physics is similar if in both cases the contact 
doping is such that the Fermi level is sufficiently be- 
low the edge of the conduction band in the sample. In 
such a case, only the tail of the distribution function is 
injected, which leads to the nondegeneracy of the elec- 
tron gas in the ballistic part of the diode. The theory 
is therefore applicable ±q quantum heterostructures with 
over-barrier transport ,£3 where current is determined by 
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a tail in the distribution function (ballistic-injection, real- 
space-transfer devices, etc.), as well as for the homodiode 
with a nondegenerate electron gas in the contacts. 

In order to simplify the problem, we assume that un- 
der the range of biases of interest, due to the large differ- 
ence in the carrier density between the contacts and the 
sample, and hence in the corresponding Debye screening 
lengths, all the band bending occurs in the ballistic base, 
and therefore the relative position of the conduction band 
and the Fermi level Ec — sp does not change in the con- 
tacts. For such a modeling, all of the potential drop takes 
place exclusively inside the ballistic base between the po- 
sitions X = and X = dva Fig. ^ and the contacts may 
be excluded from the consideration. This assumption is 
better fulfilled for the case of the hcterojunctions because 
of much higher electron densities in the contacts. 

The carriers inside the contacts are assumed to re- 
main at thermal equilibrium, and their injected part is 
distributed over the energy according to the Maxwell- 
Boltzmann distribution function at lattice temperature 
T. For the ballistic part of the diode, we suppose 

A^<rf<Ap, (1) 

with the electron wavelength and Ap the mean free 
path, so that electrons may be considered as classi- 
cal particles moving ballistically between the contacts 
and interacting with each other electrostatically. This 
regime is accessible in modern device fabrication tech- 
nologies for which the mean free path Ap may be as high 
as 10^-10^ nm in modulation-doped structures (fot-in- 
stance, in GaAs/Ala;Gai_a;As at low temperaturesEac3) 
and ~ 10'^ nm in the purest bulk material, whereas the 
Fermi wavelength is about 40 nm. 

Next we assume that the transversal size of the diode is 
sufficiently thick (much larger than the screening length 
Ld)- This allows us to treat the steady-state electro- 
static problem as a one-dimensional one in the plane 
geometry. However, to use the same one-dimensional 
consideration for the fluctuation problem, we need an 
additional justification. The fluctuating current is de- 
termined by a random transmission of discrete electron 
charges of the amount of q. Essentially, this discreteness 
of charge transmission together with randomness leads 
to the shot noise. In principle, each single electron while 
transmitted between the contacts disturbs the electric 
field and thereby interacts with other electrons of the cur- 
rent flow in both longitudinal and transversal directions. 
The electrostatic screening in such a problem is three- 
dimensional. Nevertheless, we shall treat the problem as 
a one-dimensional one considered in the plane geometry 
by averaging the fluctuations over the transversal direc- 
tions. This is justified if the average distance between 
the excess (fluctuating) carriers in transversal direction 
is much smaller than the characteristic scale of the elec- 
trostatic potential variation in that direction. This con- 
dition may be written as 

Lidvwy»i, (2) 



where is the transverse characteristic scale, n is the 
typical electron density in the ballistic region, and 5n 
its fluctuation. To estimate the order of magnitude of 
the fluctuation 5n, we use Poissonian statistics, leading 
to the relation {Sn?) ~ n/{L\d). Thus, the condition 
@) becomes nL\d ^ 1. The scale L±^ depends on the 
ratio between the longitudinal dimension d of the sam- 
ple and the Debye screening length L]j m the active re- 
gion. Fo r nondegener ate electrons, the latter is deflned as 
Ld = \J KksT / {q^n), with k being the dielectric permit- 
tivity and ks the Boltzmann constant. To estimate the 
magnitude of Lj^, we distinguish two differentpCases: (i) 
Weak screening, d Ld- For this case Lj_ ~ dj^ and the 
condition (||) becomes n ^ d^^, which for d ~ 300 nm 
requires n ^ 10^'' cm^^. (ii) Strong screening, d ^ Ld- 
For this case L± ^ Ld, and condition (0) becomes 
n '3> L^d^^. After the substitution of the expression 
for the screening length, it is seen that this condition 
becomes independent of n, although it requires a sufR- 
ciently long sample, d » q^/^nkBT) ~ 2ao{Eo/kBT). 
where oq = nh^ / {mq^) is the effective Bohr radius and 
Eo — ^^/(SKao) is the effective Rydberg energy in the 
material. For GaAs, uq « 10 nm, i?o ~ 5 meV, which 
corresponds to the temperature of about 60 K. Then for 
T ^ 10 K, d 120 nm, which is supposed to be ful- 
flUed. On another hand, the condition of strong screen- 
ing requires d^ Ld, which leads to the condition on the 
electron density 



For the same set of parameters, one gets n ^ 2 x 
10^"* cm"'^. Therefore, for both cases of weak and strong 
screening, there is a requirement on the minimal elec- 
tron density or, equivalently, on the minimal density of 
the injection current in order to use the one-dimensional 
electrostatic screening picture for the fluctuations. Oth- 
erwise, each carrier perturbs the electrostatic poten- 
tial independently and the three-dimensional approach 
is needed. On the other hand, the assumption of the 
nondegenerate electron gas restricts our approach by a 
maximum electron concentration dependent on T. For 
temperatures in the range 10-77 K, these maximal con- 
centrations are estimated to be in the range 3 x 10""^^ to 
6 X 10^^ cm^'^. These estimates show that the approach 
undertaken below covers a wide range of typical diode 
parameters: electron concentrations, diode lengths, and 
temperatures. 



III. BASIC EQUATIONS 

A semiclassical space-charge-limited transport in a bal- 
listic conductor is completely described by the electron 
distribution function F(X, Vx, t) and the electrostatic po- 
tential ip{X,t). Here, Vx is the X component of the elec- 
tron velocity and t is the time. The potential '^{X, t) 
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inside the sample is determined by the distribution of 
space charge from the Poisson equation 



with the boundary conditions 

(p{0,t) = ipL, f(d,t) = 



(4) 



(5) 



The voltage bias between the contacts U — fn. ~ fL 
is assumed to be fixed by a low-impedance external cir- 
cuit. The electron density N{X,t) at any plane X is 
determined by integrating the local electron distribution 
function over velocities 



N{X,t) 



F{X, Vx,t)dvx 



(6) 



whereas the current in the external lead is given b 



qA 

d 



d r 



Vj:F{X, Vx,t)dv^ 



dX + Co 



dU 
dt ' 



(7) 



where Cq = nA/d is a capacitance and A the cross- 
sectional area. Due to a fixed-applied-voltage condition, 
in what follows we shall neglect the last term in Eq. (|^) 
coming from the displacement current contribution. In 
addition, for simplicity, we shall omit the minus sign for 
the current, which is opposed to the direction of electron 
flow. Moreover, as will be shown below, the current is 
conserved along the sample due to the conservation of 
electron energy under ballistic motion (this is true for 
both the stationary current and its fluctuation). There- 
fore, the integration over X becomes trivial and it will 
be disregarded. 

Under ballistic motion the distribution function 
F{X,Vx,t) obeys the coUisionless kinetic equation 



OF 
'dt 



OF 

'dX 



q dip OF 
m dX dvr 



= 0, 



(8) 



where m stands for the electron effective mass. The dis- 
tribution functions of injected from the contacts electrons 
are assumed to be given as 



F{0,Vx,t)\y^^Q = FL{vx,t), 
F{d,v^,t)\y^<Q = FR{vx,t). 



(9) 



The kinetic equation (||) with the electrostatic potential 
determined self-consistently from Eqs. (^L-and (^) are 
known as the Vlasov system of equationscj describiag 
the dynamical screening of the interaction in plasma.c3 
Equation (§) may also be expressed as 



dF 
'dt 



0, 



(10) 



since F is constant along an electron trajectory, i.e., the 
distribution function at any plane X can be expressed 
through the functions Fk{vx,t), k = L,R defined at the 
boundaries. Each of these functions is considered to con- 
sist of two terms, a stationary part describing the station- 
ary injection and a time-varying stochastic component. 
Explicitly, 

Fk{vx,t)^Fk{vx) + SFk{vx,t), k^L,R. (11) 

Under nondegenerate and equilibrium conditions in the 
contacts, we assume for the stationary part of the injec- 
tion function the half-Maxwellian distribution 



Fk{vx) 



2N, 







VoV^ 



2 I 2 



(12) 



with z;^; > for k — L and < for k = R. Here, 
Nq is the density of electrons injected from the contacts 
and vo = \j2kBT jra is the thermal velocity. The contact 
distribution functions ( [l^ ) are normalized in such a way 
that the integration over a half-velocity space yields the 
density of electrons injected from the contact 



FL{vx)dvx 



v^>0 



FB,{Vx)dVx 



(13) 



l'x<0 



The stochastic terms 5Fk, k = L, R in Eq. ( |ll| ) are 
the only sources of noise under the ballistic transport 
considered here, since the electron motion between the 
contacts is noiseless. Their equal-tiipe correlation, due 
to equilibrium conditions, is given byE3 



{6Fk{vx,t)6Fk'{v'^,t)) 

= CFivx)[l-Fivx)]Skk'S{vx 



(14) 



where the constant C is determined from the normaliza- 
tion condition. Since the injected electron gas is nonde- 
generate, F <^1^ and the factor \ — F will be ignored. 

As a consequence of the fluctuations inside the con- 
tacts (whose origin is ultimately the carrier scattering 
processes), both the electron distribution function and 
electrostatic potential in the ballistic sample fluctuate, 
leading to the current fluctuations. These quantities 
will be presented as a sum of stationary and fluctuat- 
ing contributions: F[X, Vx,t) = F{X, Vx) + SF{X, Vx, t), 
N{X, t) = N{X) + SN{X, t), ip{X, t) = (p{X) + 5ip{X, t), 
and I{t) ^I + SI(t). 

Introducing the Fourier transform for the fluctuations 
of the distribution function SFt^{X, Vx) and the potential 
dipi^{X), the kinetic equation takes on the form 



iujSF^ 



dSF^ ^ q dip d&F^ ^ q dF dSip^^ 
dX m dX dvx m dvx dX 



trajectory 



with the boundary conditions at the contacts 

SF^{0,Vx)U^>^=5Ft{vx), 
SFUL,Vx)\ SF^ivx), 



0, 



(15) 



(16) 
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where 5F^ and SF'^ are the Fourier transforms of the 
stochastic functions from Eq. (|ll|). The equation for 
the fluctuating potential Sip^: is trivially obtained from 
Eqs. (|) and (|), 



(22) 



(17) 



the boundary conditions for which follows from Eq. (||), 
S^m=0, 5^1,{d)^Q. (18) 

Below we restrict ourselves to the calculation of the 
low-frequency plateau of the noise spectrum; thus one 
can omit the term proportional to uj in Eq. (|lq). It can 
be shown that this approximation is valid if the short- 
est fluctuation period in 5Fk [t) is considered to be suffi- 
ciently greater than the average electron transit time tt 
across the diode, i.e., u <C t^^. Thus, the above self- 
consistent equations completely describe the stationary 
transport and low-frequency fluctuations in the ballistic 
sample, and below we shall omit the index uj. 

It is advantageous to rescale all the variables as follows: 



__ Vx_ 
Vq ' 
TV 



X 



2iVo' 



5f = 5F 



2No' 



(19) 



In such units the basic equations contain only two di- 
mensionless parameters: (i) the length of the sample 
(or the screening parameter) A = d/L^, where = 
y^lJi37y(2g^]Vo) is the Debye screening length corre- 
sponding to the electron density 2A'o, and (ii) the applied 
voltage bias V = qU / {ksT). Below we use the dimen- 
sionless variables in all the equations. 



IV. STEADY-STATE PROBLEM 

The calculation of fluctuations in the ballistic conduc- 
tor requires the knowledge of the stationary distribution 
of electrostatic field, which, in turn, can be determined by 
solving the full steady-state problem. The self-consistent 
steady-state problem can be solved as follows. First, we 
solve the stationary coUisionless kinetic equation for the 
distribution function f{x,w) 



dx 2 dx dw 



(20) 



at a given electrostatic potential Integrating 
/(x, w) over w, we then find the electron density pro- 
file n{ip) in terms of the potential i^{x). Then we should 
solve the Poisson equation 



dx'^ 



with the boundary conditions 



(21) 



Here, we set the zero value of the potential at the left 
contact. 



A. Stationary distribution function 

To solve the stationary kinetic equation (|^, we have 
to specify the boundary conditions for this equation at a 
given ip{x). Generally, the nonstationary kinetic equa- 
tion (g) and the distribution functions (^) of injected 
electrons completely determine the nonstationary solu- 
tion f(x,w,t). However, the steady-state Equation ( [20| ) 
requires a specification of the boundary conditions for 
the distribution function of all the electrons: those in- 
jected from the contacts into the sample and those leav- 
ing the sample. Let the space charge in the sample be 
such that a potential minimum V'm occurs at x = Xm, 
which acts as a potential barrier for electrons. We define 
the total electron energy St = w'^ — ^p{x). For a given 
potential, the distribution function should consist of the 
terms originating from two electron streams injected by 
the left and right contacts. Electrons injected from each 
of the contacts fall into two groups depending on their 
injecting energies. If the initial energy is higher than the 
height of the barrier, electrons obviously reach the oppo- 
site contact and contribute to the electric current. These 
electrons are not refiected back. Note that the height 
of the barrier is different for the electrons injected from 
the left and right contacts. For those injected from the 
left, it is ipL ~ V'm = Km which is the potential mini- 
mum depth, while for those injected from the right, it is 
— V'm = V -\-Vm (see Fig. |l|). Accordingly, the lower 
bounds for the velocities of the transmitted electrons are 
given by 



wr 



= y^lpR- 1p7n = \/V„r + V. 



(23) 



Electrons of the second group, which we shall call the re- 
flected electrons, are reflected by the barrier and do not 
contribute to the current (however, both groups affect 
the electrostatic potential) . An electron from the second 
group being injected with a velocity w returns to the 
contact with the opposite velocity of the same value —w. 
Taking into account the above consideration, the elec- 
tron distribution function /(x, w) at any plane x may be 
written as 



/ = Jlj + fL,r + fR,t + fR,r 



(24) 



where the indices L and R refer to the left and right 
contacts, and the indices t and r distinguish the trans- 
mitted and reflected groups of carriers, respectively. The 
boundary conditions for these functions read 
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/L,t(0, Wc) = /l(wc) 0{wc - Wl), 
/L,r(0,Wc) = fL{Wc)9{w\ ~wl), 
/i?,t(A, Wc) = /i?,(^c) 6'(-Wc - Wr), 

fR,r{>^, Wc) = fniwc) 9{w\ - wl), 



(25) 



where Wc is the x velocity component of injected electrons 
at the contacts, 6 is the Heaviside step function, and the 
distribution function of injected electrons is determined 
by Eq. (|l2[), which in dimensionless units reads 



fLiwc) = faiwc) = -^e~ 



(26) 



We can solve now the collisionless kinetic equation 
( pO| ) explicitly for a given potential profile ip{x). In- 
deed, one can easily see that its solution is an arbi- 
trary function dependent on the total electron energy 
JF(ej) = T(w'^ — iI;{x)). The boundary conditions ( p5| ) de- 
termine the shape of this function. By using the electron- 
energy conservation law 



(27) 



where Wc and Tpk are the parameters at the contacts, we 
exclude Wc in the boundary conditions (^5|) and obtain 
the contributions in the distribution function as 



fL,tix,w) 
fR,t{x,w) 

fk,r{x,w) 



—=9(w — w, 
\/tt 



w - 



ix))e 
'w^,{x)) e 



-w' +Tp(x) — TpR 



1 



—=e{iijl{x) - e 



-w +ip{x)—ipk 



(28a) 
(28b) 
(28c) 



Xm and the lower sign for the right side of the potential 
minimum Xm < x < X. It is seen, that the obtained 
distribution function is discontinuous on w at the points 
where w = w,(a;) (see also Fig. || discussed below). It 
is not surprising, since only a discontinuous solution can 
satisfy the first-order equation (^o|) and simultaneously 
two different arbitrary functions given at the boundaries. 
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FIG. 2. Typical electron trajectories in the phase space 
{x, w) for different electron groups: L and R refer to the car- 
riers originated from the left and right contacts, and t and r 
refer to the transmitted and reflected groups of carriers. The 
separating curves are the critical velocities ±ui*(a;), which in- 
tersect at the point of the potential minimum {xm,0). The 
results are for A=30, V=10. 



where k — L, R, and the functions fL.t and /rj for 
the transmitted electrons are defined in the whole range 
< X < A, whereas the expressions for the reflected elec- 
trons fL,r and ffi^r are valid in the intervals < x < Xm 
and Xm < X < X, respectively. In Eqs. ( ^ we have 
introduced the quantity 



W^:{x) = \/i>{x) - Ipn 



(29) 



which has a meaning of the maximal velocity of reflected 
electrons at a point x. For the sake of clarity, in Fig. ^ we 
show the electron trajectories in the phase space (x, w) 
corresponding to different electron groups. It is worth 
to stress that the distributions (|2|) depend on the local 
potential tpi^) and the potential minimum ipm as well, 
i.e., the distribution function depends non locally on the 
potential profile. 

Summing up all the contributions (p8|), the total dis- 
tribution function takes on the form 



1 



fix,w) ^ —=t 



— IpR 



w > =Fw*(x), 
w < =Fu;*(a;). 



(30) 



Here, and throughout the paper, we shall use the upper 
sign for the left side of the potential minimum < x < 



B. Electron density 

The distribution function ( |30| ) allows us to find the 
electron density at a slice x as 



n(x) = —;=e 
a/tt 



ijj^x) 



{x) 



-4'R 



dw 



^10^ (x) 



dw 



= ie'A(^)(e-'^^{l ± erfK(a;)]} 

+e-'^«{l^erfK(a;)]}), (31) 

where erf(a;) — (2/^/??) /o^e~" du stands for the error 
function. By using the values for the potential (|2^) at 
the contacts and denoting 



(32) 



the electron density can be written as a function of ip, 



6 



(33) 



where, as before, the upper sign apphes in the interval 
< a; < Xm and the lower sign applies in the interval 
Xm < X < X. Note that in equilibrium, V = 0, Pi = 2, 
(32 = 0, the Boltzmann distribution n{x) = e'^'^-' is re- 
covered throughout the sample. Furthermore, Eq. ( |3^ ) 
is valid for single- injection (vacuum) diode, assuming 

In the following we shall use the shifted potential mea- 
sured from the minimum 

r]{x) = ipix) - Iprn, (34) 

and Eq. (|3^ ) in terms of the new variable r] becomes 

n{f]) ^ n^ns"^ [1± Peiiy/f]\, (35) 

where — ^Pie"^'" is the electron density at the po- 
tential minimum, and 



/3^|=tanh(| 



C. Steady-state electrostatic potential 



(36) 



a; = < 



dr] 



< X < Xm 

Xjyi ^ X \ 



(41) 



where the boimdary conditions for ri{x) are 

T](Q) = riL = Vm, i-i{\)=m = y-m + V. 



(42) 



For the given V^,A, the only unknown parameter in Eqs. 
( ^ ) is the potential minimum Vm- The latter is found 
by matching Eq. ( p| ) at x = Xm, where r]{xm) — 0, and 
one gets 



where 



dr] 



. Ihy{r]) 



V^(V) + V 



dr] 



Xn 



XV^r 



(43) 



(44) 



is the screening parameter renormalized to the electron 
density at the potential minimum rather than to the con- 
tact electron density as before. 



Having found the analytical expression for n{ip), we 
have to use it to solve the Poisson equation (|l]). Multi- 
plying both sides of Eq. (|2^) by dtp/dx and integrating, 
one gets 



dx 



n{ip)dip, 



(37) 



where we have used the property of the potential min- 
imum (dip / dx)\x='x^ = 0. Changing to the shifted- 
potential variable r] and carrying out the integration, one 
gets 



where l/l"^ 



Pie' 



and the function 



(38) 



(39) 



depends on the applied voltage V through /3. Taking into 
account drj/dx < for < a; < Xm and drj/dx > for 
Xm < X < X, the electric field is given by 



h^ir]) - e" - 1 ± /3 e^crfV^ - 



dx 



\J hyir])/e,n, < X < Xm 
-\J hvi'ri)/^m, Xm < X < X 



(40) 



which is measured in units of kgT/qL^. Integrating Eq. 
(^), one obtains the distribution of the potential in an 
implicit form. 



D. Steady-state current 

This brief description of the steady state is then com- 
pleted by the expression for the stationary current. Sub- 
stituting the distribution function into Eq. (M) and chang- 
ing the variables with the help of Eq. l\2T\ ) as wdw = 
Wcdwc, one obtains 



fL,t{x, w)wdw 



-w* {x) 



2^/^Ic 



fR,^t{x,w)wdw 

fL.t{0, Wc) WcdWc 

fR,tiX,Wc)Wcdw. 



where 



-.qNoVoA = qNovA 



(45) 



(46) 



is the emission current from each contact (limiting value 
for the total current at V oo, Vm 0), and v — 
vq/^/tt ~ y/2kBT/ (TTm) is the average velocity of the 
injected electrons with the half-Maxwellian distribution. 
Only the part of the distribution function corresponding 
to the transmitted electrons has been taken into account, 



7 



since the reflected carriers gives no contribution to the 
current. [This is in contrast to the case of the calculation 
of the electron density ( |3l| ) for which both transmitted 
and reflected carriers contribute.] It is seen from Eq. 
( ^ ) that the current is the same for any section x of 
the sample, given by its value at the injected contacts. 
Substituting the functions ( p5| ) into Eq. (^) and carry- 
ing out the integration, we obtain the current as a sum 
of two opposing currents: Ilr and Irl caused by the 
injection from the left and right contacts, respectively. 



/ = IrC- 



LR 



I 



RL- 



(47) 



The formula for the current may be written through the 
electron density at the potential minimum, that is. 



/ — 2n.„iIcP — qN„ivA tanh 



2k„T 



(48) 



where Nm='2Nonm- This formula justifies the usage of 
the term "virtual cathode" referred to the location of the 
potential minimum, since it is seen that the current is 
determined by the injection of the electron density Nm 
from the virtual cathode. The additional tanh( ) factor 
takes into account the injection in the opposite direction, 
and it tends to 1 a,t qll ^ ksT. (For the vacuum diode 
case, this factor is set to 1 because of only one injecting 
contact.) 

Summarizing this section, we note that the above re- 
lations solve completely the steady-state problem for the 
ballistic two-terminal conductor: Eqs. (|41 ) determine the 
distribution of the potential across the diode in an im- 
plicit form, and Eqs. (^3|) and (^) determine the current- 
voltage characteristics. Note that in Eq. (|4^) the cur- 
rent depends on voltage through both the explicit term 
and the potential minimum Vrm which is a func- 
tion of voltage. Equations (|33|), (|9|)-(|3|), (|^ may 
be viewed as an extension of the EcypUangmuir theory 
for a single-injection vacuum diodecj l3 to the double- 
injection case. The Fry-Langmuir formulas are obtained 
by setting (32 = 1, Irl=0. 



E. Results 

Figure || shows the typical spatial distributions of the 
potential ■0, electric field E, and electron density n along 
the diode obtained from Eqs. (||), (||), and (H). 

With the aim to compare ourJieory with the results 
of the Monte Carlo simulations Jl3 we present the spatial 
profiles for the value of A=30.9 and various applied bi- 
ases V. As it is seen from the figure, the agreement is 
excellent for all the quantities. 

The space-charge-limited conduction is characterized 
by a strong transport inhomogeneity in the ballistic re- 
gion and by the presence of the potential minimum [Fig. 
^(a)] due to the injected space charge. The minimum acts 
as a barrier for the electrons moving in both directions. 



Its magnitude progressively decreases as the applied bias 
is increased, simultaneously shifting towards the left con- 
tact on which the potential is fixed. 




FIG. 3. Spatial profiles for the normalized quantities: (a) 
potential ip, (b) electric field E, and (c) electron density n 
(all solid lines) for A=30.9 and several applied biases V . The 
corresponding units are ksT/q, ksT/qL^, and 2A^o. The re- 
sults are shown to be in excellent agreement with the Monte 
Carlo simulations (Ref. 48) (symbols). 

The obtained solutions are determined by two dimen- 
sionless parameters: A and V. The spatial distributions, 
however, may be presented in a more universal form by 
using for scaling the potential minimum parameters. We 
define the new coordinate X ^ ^ Xm)l(-m, where the 
characteristic length Ira = (2nm)^^^^, dependent on the 
electron density at the potential minimum, has been in- 
troduced in Sec. tVC . This is equivalent to scale the 
original coordinate X in units of the screening length re- 
ferred to the electron density at the potential minimum 
rather than to the contact electron density. In such a 
unit, the parameter A is scaled away from the equation 
for the potential, remaining only in the upper and lower 
bounds of the function variation. Explicitly, Eq. (|4T 
becomes 



X = 



- Io[dv/\/hv{v)], -Xrnl^rn < X < 0, 

< X < (A - x™)/£„ 



(49) 
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FIG. 5. Current and its components coming from two op- 
posite electron flows I = Ilr — Irl (in units of Ic) vs bias 
V for A=50. The height of the potential barrier Vm, and its 
location XmjX are also shown. 



FIG. 4. Profiles for the steady-state quantities in units of 
the potential-minimum parameters: (a) potential, (b) electric 
field, and (c) electron density. The limiting universal profiles 
at y ^ oo are shown by solid lines. 

Therefore, all the solutions may be presented as a one- 
parameter family of curves dependent on the applied bias 
V only. Moreover, at high-voltage limit V > 5, /3 ^ 1, 
the functions hy{r]) become independent of bias, and the 
spatial distributions tend to the limiting universal pro- 
files for each quantity which are free from any param- 
eter. This is valid for all the spatial characteristics as 
it is seen from Fig. ^ where the potential rj = ip — 
the electric field Eim, and the electron density n/rim are 
plotted. Moreover, at this limit the part of each pro- 
file at X < tends to vanish (the potential minimum 
approaches the left contact), which leads to the validity 
of the virtual-cathode approximation with the bound- 
ary condition i?(0)=0. The universal profiles obey the 
asymptotic behavior at x ^ oo: 77 (x) = fax^^'^j Elm = 
-ax^/^ n/um = |ax"^/^ where a = (3/7r)i/3 w 0.9847. 
Going back from x to the x coordinate, and using the 
Child law (^4|), which will be discussed below, one can 
obtain the asymptotic formula for the potential profile 
rjix) = V{x/ X)'^^^ , X 00. The latter is valid not 
only for the present nondegenerate-electron-gas model, 
but for an arbitrary distribution function of the injected 
electrons, provided V —^ 00, X —^ 00 (virtual-cathode 
approximation) .EJ The related formulas for E{x), n{x) 
may also be obtained by taking the derivatives. 



The choice of the potential minimum parameters as ref- 
erence coordinates is of traditional use in vacuum-diode 
literatureBQE^ Since only one contact (cathode) is con- 
sidered as injected for these diodes, /3=1 for any bias, 
and the universal potential profile independent of the 
diode parameters is obtained for any bias, as it was tabu- 
lated in the original work by Langmuir.c£l For the case of 
the two-terminal semiconductor diode, that universality 
is broken at low and moderate biases due to the contri- 
bution to the current from the second injecting contact, 
but it is recovered however at high biases V —^ 00 when 
the influence of the second contact becomes negligible, 
as it is demonstrated in Fig. We remark addition- 
ally that the virtual-cathode approximation is only valid 
when besides V ^ 00 another condition is fulfilled si- 
multaneously, V < Vcr- Otherwise, the transport is no 
longer limited by the space charge, the current saturates 
at / = Jc, and the value of the electric field at the left in- 
jecting contact is no longer zero, £'(0) < 0. This change 
in the transport regime is clearly seen in Fig. ^ where 
the current and its components coming from two opposite 
electron flows / = Ilr — Irl versus bias V are plotted 
for a particular value of A. It is seen that the current is 
an increasing function of the bias up to the critical value 
Vcr, after which it is saturated at I ~ Ic- At that point 
the potential minimum vanishes. For V ^ 5 the contri- 
bution to the current from the right-contact electrons is 
also essential. 
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FIG. 6. Current-voltage characteristics for different lev- 
els of screening A obtained as solutions of the present theory 
(diamonds). For comparison, the approximate solutions are 
shown: linear dependences given by Eq. @ (solid); Lang- 
muir formula given by Eq. ( |54[ ) (dots); the Child 3/2-power 
law given by Eq. ( js^ ) (dashes); parameter /3, shown by long 
dashes, when approaching to 1 indicates the bias {V ~ 5) 
over which the effect of the second contact on I-V curves be- 
comes negligible. Monte Carlo simulation results (Ref. 48) for 
A=30.9 are shown by squares. 

The I-V curves for different levels of screening are 
shown in Fig. We have checked that the obtained so- 
lutions are in excellent agreement with the Monte Carlo 
simulationsca (the case of A=:30.9 is compared in the fig- 
ure). 

The analysis shows the following behavior. At low bi- 
ases, the I-V curves are linear for all A despite the fact 
that the transport is space-charge-limited. The curves 
for this case are described by 



'lin 



IcVe- 



V<1, 



(50) 



where V^ is the equilibrium value of the potential mini- 
mum whose value depends on A. In the range 1 ^ ^ 10 
the I-V curves deviate to sublinear dependence. At high 
biases, starting approximately at y ^ 5 where /9 ^ 1, 
the effect of injection from the second contact becomes 
negligible, Irl ^ Ilr (see also Fig. ||). Furthermore, for 
V 3> Vm ^ 1 the analytical solution may be found. In 
this regime, the function (jf) may be approximated by 
leading-order terms of a series expansion in a similar way 
as in the case of a vacuum diodecZl 

/i+(7?) « 2V?7A- 1, ry-^oo. (51) 

In this regime from Eq. ( pl| ) one can write 

dr] 



A 



V4 (2 VW^- 1)1/2' 



(52) 



from wkiph|-by using £rn ~ 1/VJ follows the Langmuir 
formulaB'E§EZl 



/l 



ang 



8 {V + Kn)3/2 



(53) 



In Fig. g we present the curves calculated from this for- 
mula, and they are seen to describe accurately the I-V 
characteristics for the highest biases. For higher A, the 
range of biases where this formula may apply is wider. In 
the asymptotic limit ^ oo, A — > oo, one may neglect 
Xjn and Vm as compared to A and V, respectively, and 
one obtains the Child 3/2-power law, which is free from 
the potential minimum parameters 



^Child = g"^-^' 



1/3/2 



(54) 



It is seen from the figure, that this asymptotic formula 
accurately describes the I-V curves only at very high val- 
ues of the parameters: A 10^, V ^ 10'^. However, as 
we have discussed earlier, there is a relevant difference 
between the semiconductor and vacuum ballistic diodes. 
In vacuum diodes the applied voltage may be quite large 
without breaking down the ballistic transport regime. In 
contrast, in solids, electrons even for a pure material in- 
teract with a lattice. Under a low-bias regime this in- 
teraction is weak, but it becomes quite strong at high 
biases due to the significant increase of the electron en- 
ergy. For instance, the threshold for the optical phonon 
generation in GaAs is about 0.036 eV, which corresponds 
to « 40 at r ~ lOK. Thus, one cannot bias the sam- 
ple to the voltage more that that value, since a strong 
interaction with the lattice will break down the ballistic 
regime. The allowed range of biases is typically restricted 
hy U ^ BOksT/q. Then, for real structures the ballistic 
length A is well below 100. Therefore, the Child 3/2- 
power law is hard to achieve in semiconductor ballistic 
n-i-n diodes, and one should use the full set of formu- 
las described in the present paper from which follows the 
linear or sublinear I-V dependences in a wide range of 
biases even under a strong limitation of transport by a 
space charge. 

Finally, Fig. |^ illustrates the stationary electron distri- 
bution function over velocities f{w) at different sections 
of the diode for several biases V. The distribution func- 
tions are discontinuous ai w — 'w*{x)sgii{x — Xm)i as 



discussed in Sec. |VA|. It is interesting to note that at 
high biases the arriving electrons at the right (receiv- 
ing) contact exhibit a sharp peak separated in energy 
from the intrinsic contact electrons [f^ig R(d)] and thus 
may be distinguished in an experiment While the in- 
jected carriers are uncorrelated, electrons arriving at the 
receiving contact that belong to that peak exhibit corre- 
lations in en ergy. This interesting result will be discussed 
in Sec. IVH. 
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Shie) = V^IcSMe), k = L,R, (58) 

where SIk is the low-frequency Fourier component of the 
injection current fluctuation, and Ic is the mean emis- 
sion current defined by Eq. (^6|). The correlator for SIk 
is obtained from that for 6Fk given by Eq. (|lj) , and one 
gets 

{Slkie)dlk'{s')) = 2qI,Afe-'dkk'6ie - e'), (59) 

with A/ the frequency bandwidth. The obtained corre- 
lator shows that the electrons with different energies are 
uncorrelated, which is a consequence of the Poissonian 
injection statistics. The fluctuations at the left and right 
contacts are assumed to be uncorrelated as well. 



FIG. 7. Electron distribution function over dimensionless 
velocities /(w) at different sections x of the diode for several 

biases V and A=30. A. Injected electron-density fluctuations 



V. FLUCTUATION PROBLEM 

We will find here the fluctuations of the distribution 
function, electron density, electrostatic potential and cur- 
rent in the ballistic region of the diode, which are caused 
by the fluctuations in the contacts. To solve the fluctu- 
ation problem (p^)-(|l8|) self-consistently, we undertake 
the same approach as used above for the steady-state 
problem. First, we calculate the fluctuation of the dis- 
tribution function (5/ in a given electrostatic potential 
'0(x) -|- S'ip{x) by solving the perturbed kinetic equation 



dSf 1 di} 85 f 1 df d5^ 
dx 2 dx dw 2 dw dx 



= 0. 



(55) 



The solutions for 5f{x, w) for different groups of elec- 
trons are derived in the Appendix, where we also present 
the corresponding electron density fluctuations Sn{x) ob- 
tained by integration over velocities. The fluctuations 
6n, which are the functions of 6ijj and the contact fluctu- 
ations (5/l and Sfu, should then be substituted into the 
perturbed Poisson equation 



dx'^ 



= Sn{x) 



with the boundary conditions 

(5V'(0) = 6ip{X) = 



(56) 



(57) 



to find the self-consistent fluctuations of the potential 5ip. 

The fluctuations of the distribution functions of in- 
jected electrons Sf^, SJr are supposed to be given by 
the correlator (p^. It is advantageous, however, to ex- 
press them through the injected current fluctuations. For 



each injection energy e 
by 



their relationship is given 



The electron-density fluctuations at a slice x caused by 
the stochastic injection from the contacts is obtained by 
summing up all the contributions (A15) derived in the 
Appendix. In terms of the injected current fluctuations 
(pq), we obtain the following expression: 



(5n"J(x) 



5 hie) de 



( r^^-^r. Sh{e)de 

< X < Xn 



(60) 



Xm < X < A. 



B. Induced electron-density fluctuations 



The electron-density fluctuations induced at a slice x 
by the fluctuations of the pote ntial is obtained by sum- 
ming up all the contributions ( Al6| ), and one gets 



6n''^'^{x)^n{x)5i){x)± 



J 



2y/TTW^,{x 



■[(5V(x)-<5V™], (61) 



where J = I / Ic, and the upper sign applies in the inter- 
val < a; < Xm and the lower sign applies in the interval 
Xm < X < X. This term along with the term ( |60[ ) should 
then be used in the perturbed Poisson equation. 
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C. Current fluctuations 

The expression for the fluctuation of the current in any 
section of the sample is given by 



/oo 
-oo 



w) w dw. 



(62) 



Now we have to substitute into Eq. (|6^ ) the fluctuation of 
the distribution function, which is convenient to consider 
here as a sum of the homogeneous and nonhomogencous 
parts of the solution of the kinetic equation [see Eq. ( A4) 
in the Appendix]. The contribution of the nonhomoge- 
neous term is zero, whic h can be easily checked by direct 
integration of Eq. (A7). The homogeneous term con- 
sists of the transmitted and reflected parts given by Eqs. 
(A6). Again, the reflected electrons give zero contribu- 
tion to the current fluctuations, since the functions ( A6b) 
are even on w, so that the integrand ( |62[ ) is an odd func- 
tion and its integration from — oo to oo yields zero. The 
only nonzero contribution comes from the terms (A6a) 
for transmitted perturbing electrons. Substituting them 
into Eq. ( |6^ ) and changing the variable of integration 
from w io Wc, we obtain 



51 = 2v^/, 



C30 

SfLiWc)WcdWc 

+ I 5fR{wc)wcdwc 



+Ic e"^™ (e-"^^ 



-4'R 



(63) 



from which it is seen that the current fluctuation is inde- 
pendent of the position x. By using the definition for the 
injected current fluctuation (BSl) and the formula (47) for 



the average current, the final expression for the current 
fluctuation takes on the form 



SI ^ / 5lL{e)de- 



SlRie)de-ISVm, (64) 



where SVm = —S'>p{xm) is the potential minimum fluc- 
tuation. 61 depends on the magnitude of the fluctuat- 
ing potential barrier irrespective of its random location. 
This is a consequence of the current conservation along 
the diode. 



Equation (|64) is a central one, which determines the 
fluctuation of the transmitted current through the fluc- 
tuations injected from the contacts. The first two terms 
in the rhs represent the current fluctuations transmit- 
ted directly to the opposite contact from the left and 
right contacts, respectively. Since the injected electrons 
of different energies are uncorrelated, they give the full 
shot noise. It is the last term —ISVm, caused by the 
self-consistent potential fluctuation (long-range Coulomb 
correlations), that compensates the current fluctuation 
and may result in the noise reduction. We note, first, 
that it is proportional to the current and thus exists only 
under nonequilibrium conditions. Second, it depends on 



the potential barrier fluctuation SVm- When the barrier 
does not appear under certain conditions, all the injected 
fluctuations are transmitted to the opposite contact and 
the noise of the transmitted current is expected to be 
the same as thet for the injected carriers, i.e., the full 
Poissonian shot noise. The compensating behavior may 
occur only when the potential barrier is present. Notice 
that the contributions of the left- and right-injected fluc- 
tuations are of the opposite sign, i.e., SIl > increases 
the fluctuation of the transmitted current, while SIr > 
decreases it. 

Among all the injecting perturbing electrons, only 
those able to pass over the potential barrier contribute to 
the transmitted current fluctuation. This fact is reflected 
in the lower integration limits that contain the height of 
the potential barrier. In contrast, all the injected elec- 
trons contribute to the potential barrier fluctuations, and 
thereby participate in the compensation effect, as it will 
be shown in the next section. 



D. Self-consistent potential fluctuations 

We flnd the potential barrier fluctuation SVm, which 
is of prime interest, from the linearized Poisson equation 
( ^ ) for the potential fluctuations Sip. By substituting 
the electron-density fluctuations Sn consisting of the in- 
jected and induced contributions found in Sees. |V A| and 



|VB| , we obtain for the self-consistent potential fluctua- 
tions 



d^Sip 
dx'^ 



= n{x) 5ip{x) ± 
-K(5n"^(a;). 



J 



2y/TTW^,{x] 



[Sip{x) - S'ip.„ 



(65) 



This is a second-order nonhomogeneous differential equa- 
tion with spatially dependent coefficients, where the term 
Sn"'-^ , dependent on the fluctuations at the contacts SIk 
[see Eq. (^0|)], plays the role of a stochastic noise source. 
To find its solution in a general form is a complicated 
problem. In addition, we remark that the term with 
l/u'*(x) is singular at the potential barrier minimum 
X = which produces an additional difficulty. Never- 
theless, we will show that it can be solved exactly with- 
out any approximation. First of all, it is advantageous 
to introduce a new stochastic quantity 



Sr]x = 5ip{x) - Sipr, 



(66) 



which is the the potential fluctuation at a slice x mea- 
sured from the fluctuating potential minimum. Thus, due 
to our choice, at the potential minimum Srjx,„—0, where 
Xm = x'^ + Sxm is a stochastic location of the poten- 
tial minimum fluctuating around its steady-state position 
Xm- The latter fluctuation, however, may be neglected. 
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since it is only of second order in respect to the poten- 
tial fluctuations, because of the property of the minimum 
tp'{x'^) — 0. Thus, one gets the stochastic differential 
equation 



J 



Silx 



(67) 



The boundary conditions for this equation follows from 
Eqs. @ and ^ 



(68) 



SrjL = SrjR = -Hra- 



Since the potential 5rjx is referenced to the fluctuating 
minimum, its values on the contacts are not zero, while 
in a stationary frame they are zero due to a fixed-applied- 
voltage conditions. 

To find the solution of Eq. (^), we use a method we 
have recently applied for a stochastic drift-diffusion equa- 
tion which has a similar form.E3 Essentially, this method 
is based on the possibility of finding two (arbitrary) lin- 
early independent solutions of the corresponding homo- 
geneous equation L5rjx — 0, which can further be used 
to construct the solution for the nonhomogeneous equa- 
tion satisfying the appropriate boundary conditions. One 
of the solutions is proportional to (dip/dx), which can 
be seen by differentiating the Poisson equation ( pTj ) and 
comparing the result with Eq. (|6^ ) with zero rhs. For 
convenience, we take it as E{x) = — {dip/dx), so the soIlu 
tion coincides with the electric field profile. In general,Lj 
the second solution can be obtained from the first one by 
using the formula u{x) = E{x) j f~,\W {y) / E"^ {yy\dy , where 
W{x) = E{x)u'{x) — E'{x)u{x) is the Wronskian, C is an 
arbitrary constant, the prime stands for the derivative, 
and E{x) ^ 0, Va; is assumed. However, this formula can- 
not be applied for our problem, since E{x)=0 precisely 
at the point of the potential minimum and the integral 
diverges. Alternatively, we use another formula for the 
second solution u which has no divergence in the whole 
region. Explicitly, 



u{x) 



W{x) 
' E'{x) 



E{x) 



c 



W{y)Q{y) 
[E'{y)Y 



dy, 



(69) 



where the function Q{x) — ~n{x) =F J/ ^4TTri{x) is a 
free term in the operator L, and the necessary condition 
E'{x) ^ is fulfilled. Next we notice that the differ- 
ential operator L given by Eq. (|67| ) does not contain the 
term with the first derivative, which leads to the constant 
Wronskian. The value of this constant is not actually im- 
portant, since it will be canceled as will be seen below so 
we take W{x) = 1. The arbitrary constant C in Eq. ( |69| ) 
does not influence on the final results. It is convenient, 
however, to define it by the conditions u{0)—u{X)—0 at 
the ends of the diode, which correspond to the homoge- 
neous boundary conditions for the Green functions of the 



operator L and provide the most compact intermediate 
expressions. To satisfy the zero boundary conditions on 
both ends of the diode, one can take the function u{x) 
as consisting of two branches. As a result, we obtain the 
following expression: 



u{x) 



1 




n{x) 

My) 



E{x) 
n{y) 



"^(y) 

^ Jv{y) - n{y) 



dy- 



1 



n^(y) 



dy 



ulEl 
1 

nnEu '' 



< X < Xr. 



< X < X, 



(70) 



where i'{x) = 1/ y^ATrr]{x) and n{x) and E{x) are the 
steady-state spatial profiles of the electron density and 
electric field, which take the values at the left and right 
contacts uljEl and nn^En^ respectively. The function 
u{x) > is continuous in the entire region < a; < A, 
including the point of the potential minimum, where it 
takes the value u{x,rL) — ^jnm- At that point, however, 
it has an infinite derivative, which is a consequence of the 
zero of the field. 

The general solution of Eq. (|67|), satisfying the 
boundary conditions ( |6^ ) and the conditions E{x,n)~^^ 
^Vxm—^j then reads 



Sr]x = E{x) / u{y)Ssy dy + u{x) 



-SljL 



E{x) 
El 



Q < X < Xr, 



E{x) j u{y)5sy dy - u{x) 
E{x 



+SriL 



Ep 



Xttt X 



E{y)Ssy dy 

(71a) 
E{y)Ssy dy 

(71b) 



where 6sx = n{x)SriL — (5n*"^(x) is the nonhomogeneous 
part of Eq. (|6^). Thus, one can find the potential fiuctua- 
tion Srjx at any section x of the sample. In particular, its 
value at the boundaries yields the potential barrier fluc- 
tuation dVrn = 6rjL. We flnd the unknown 5rjL from the 
continuity condition on the derivative dSrj/dx at a; = Xm- 



6r]i 



1 


1 ' 




E^' 




-I 



u{x)dsx dx. 



(72) 



Now recalling that drj^ has entered also in ds, we obtain 



5Vm = 



u{x)Sn^"-\x) dx, 



with 



A = 



1 

E^ 



1 

E^ 



u{x)n{x) dx. 



(73) 



(74) 



The last integral can further be reduced by substituting 
n — —dE/dx and the expression for u{x) given by Eq. 
([70|). Integrating by parts, one gets the simple formula 
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A 1 1 
A = - H . 

2 El Eji 



(75) 



The obtained analytical expression ( [73| ) with the param- 
eter A given by Eq. (j?^) yields the fluctuation of the 
barrier height in terms of the spatially distributed "noise 
source" Sn"''^ (x) caused by the random injection from 
the contacts. The weight function u{x) shows the rela- 
tive contributions of the "noise sources" to the potential 
barrier fluctuations. Its behavior is illustrated in Fig. ||, 
where we present u normalized to as a function of 

the coordinate x — ~ Xm)/im- In such a scaling for 
a fixed voltage, u{'x)nm is almost independent of A with 
a slight deviation at the ends of the function extension. 
An interesting property of those functions for different 
biases is that they cross the curve nm/n{x) (the inverse 
universal density profile as discussed in Sec. IVE) at two 
characteristic points: the potential minimum x—^ where 
du/dx — oo, and at the maximum of u{x) {see Fig. ||). 
The latter point has significance in that the electron- 
density fluctuations there have the largest influence on 
the potential-barrier fluctuations. It is worth noting that 
the maximum contribution to 6Vm does not come from 
the potential minimum location, as it would seem intu- 
itively. 




FIG. 8. Function u{x) which shows the relative contribu- 
tions of the noise sources Sn^"-' (x) to the potential barrier 
fluctuations for A=100 and for several biases V . The poten- 
tial minimum is located at x=0. 



E. Current noise spectral density 

Substituting the obtained formula (^3|) for 5Vm into 
Eq. (c4), we obtain the current fluctuation as 



^L{e)5lL{e)de + I ^R{e)SlR{e) de 

JQ 

i -2J jf^ K{x,e)dx, e<Vm, 
[l- J Jo K{x, e) dx, e >Vm, 



(76) 
(77) 



-2 J /jl K{x, e-V)dx, e< Vm + V, 



t -1 - J J(, K{x, £-V)dx, £> Vrn 



(78) 



where K{x,e) — u{x)/[2^/nA^/e + ^(J{x)], and and 
x*j^ are found from e=—'ip{x*]^)=V — ip{x*j^). The func- 
tions 7fe(e) introduced for each contact have the mean- 
ing of current fluctuation transfer functions, since they 
represent the ratio of the transmitted current fluctua- 
tion to the injected current fluctuation for a particular 
injection energy e. The terms proportional to the cur- 
rent J originate from the potential minimum fluctua- 
tions, whereas the constant contributions (±1) represent 
the direct transmission of fluctuations to the opposite 
contact. 

Equation (|7^) leads to the spectral density of current 
fluctuations 



Si - 2qlc 



de. 



(79) 



This equation with 7^(2) given by formulas (||) and (|7|) 
is the final result of our derivations. It allows us to obtain 
the current-noise spectral density, for the given level of 
screening A and applied voltage V, from the steady-state 
distributions of the potential 'ip{x), electric field E{x), 
and electron density n{x) by direct integration. Thus, 
the current-noise level is directly related to the transport 
inhomogeneity in the system. Note that the obtained 
formulas are exact for biases ranging from thermal to 
shot-noise limits under a space-charge-limited transport 
conditions. 

For practical calculations of the transfer functions 
7fc(e), one may integrate by parts the function K in for- 
mulas (frzl ) and (^, which leads to the following expres- 
sions corresponding to each group of carriers: 



7L,r(e) 



2Ar. 



■dri, e < 0, 



(80a) 



7L,t(£) = 1 - 



2A, 



[hyivW^' 



■drj 



^"-^'^ - 2A™ [/.+ (^)]3/2 

7i?,t(£) = 7L,t(e-l^)-2, e>0, 



dr]j , e > 0, (80b) 

dr], e < 0, (80c) 
(80d) 



where e — e ~ Vm is the injection electron energy refer- 
enced from the potential minimum. 



1 



1 



a/tt 



(81) 

(82) 
(83) 
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Formulas (BOa) and ( pOb| ) with /3=1 correspond to the 
formulas for a vacuum diode found by North within dif- 
ferent approach [see Eqs. (31) and (38) of Ref. ||. 



F. Nyquist equilibrium noise 

In equilibrium, / — > 0, the compensating term ISVm 
in Eq. (|6^ ) vanishes, and, comparing with Eq. ([76|), the 
transfer functions are simply the step functions with a 
step at the barrier height: ■'y'j^{e) — 6{e — Vm), 7^ (e) — 
—9{e — Vm)- This means that only electrons able to pass 
over the barrier contribute to the equilibrium (thermal) 
noise. For this case, one can easily obtain the Nyquist 
noise formula 



(84) 



where go — dl / dU\u^o is the zero-bias small-signal con- 
ductance. [To find the conductance we have made use of 
Eq. (^).] Both electron streams, from the left and right 
contacts, equally contribute to the Nyquist noise. The 
space-charge effect on the equilibrium noise is present in 
the dependence of go on the potential minimum Vm ■ 



G. Noise-reduction factor 



The obtained formula ([73|) for the current-noise spec- 
tral density Sj, which accounts for the long-range 
Coulomb correlations, may be compared with the un- 
correlated value through the so-called noise reduction 
factor. Out of equilibrium, if one neglects the term 
IdVm in Eq. ( p^ ) , which is responsible for the long-range 
Coulomb correlations between the carriers, one obtains 
7E"™"(£) = d{e-Vm), 7"(£) = -Oie-Vm-V), which 
leads to 



Sr""'^2qilLR 



Irl) = 2<7Jcoth(y/2) 
~ 2ql, V>5, 



(85) 



which is nothing more than the Poissonian noise of two 
uncorrelated streams of carriers opposite each other (at 
high voltages the contribution from the right-contact 
stream becomes negligible). It is reasonable, therefore, 
to define the noise-reduction factor by 



Si 



Si 



^unc 



2q/coth(t//2) 



(86) 



By this definition, hnth the thermal noise and shot noise 
limits are included. □ 

Figure ^ shows F versus applied voltage V for var- 
ious screening parameters A. At low values of A, the 
noise-reduction effect is weak, F w 1. As A increases, 
the noise becomes substantially reduced in the range of 
biases ksT qU < qUcr, where Ucr is a critical voltage 
for which the potential minimum vanishes (its value is 



a function of A). At U > Ucr the full shot-noise level 
is abruptly recovered. This sharp increase in the noise 
intensity when observed in an experiment would indicate 
on the disappearance of the potential barrier controlling 
the current. 
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FIG. 9. Current-noise reduction factor F vs bias U for dif- 
ferent levels of screening X—d/Lo (solid). For comparison. 
North's asymptotic solution given by Eq. is shown (dots) . 
For the case of A=30.9, the results are shown to be in excellent 
agreement with the Monte Carlo simulations (Ref. 5) (trian- 
gles). 
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FIG. 10. Contributions to the current-noise-reduction fac- 
tor F corresponding to different electron groups for the case 
of A=50. L and R refer to the left and right contacts, and 
t and r distinguish the transmitted and reflected groups of 
carriers. North's asymptotic solution is shown by dots. 

We have compared our results for the noise reduc- 
tion factor psvith those obtained by the Monte Carlo 
simulations .0 The agreement was found to be perfect 
within numerical uncertainty of the Monte Carlo algo- 
rithm, as it is seen from Fig. |^ where we show such a com- 
parison for A=30.9. The agreement for the noise charac- 
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teristics, as well as for the steady-state spatial profiles 
and I-V curves, indicates the correspondence between 
our kinetic theory and the Monte Carlo model used in 
Refs. |,| and ||. 

An advantage of our analytical approach is that, in 
addition to the net noise characteristics, one may distin- 
guish the relative contributions to the noise from different 
groups of carriers. In Fig. |l^ we present the results for 
the noise-reduction factor F as a sum of four contribu- 
tions. It is seen that in equilibrium only the transmitted 
electrons contribute to the noise (equally from the left 
and right contacts). In the range 1 ^ V ^ 10, the con- 
tribution from the reflected carriers becomes appreciable 
with a maximum at V ^ 3. At higher voltages, as the 
potential barrier progressively decreases, the role of the 
reflected carriers becomes less important. The contribu- 
tion of the right-contact transmitted electrons is negligi- 
ble at y ^ 5, as for the stationary I-V characteristics. As 
a result, in the high- voltage limit, only the left-contact 
transmitted electrons contribute to the noise. This fact 
can be taken into account in analyzing the asymptotic be- 
havior of the noise-reduction factor at high- voltage limit. 
In this limit the main contribution to the current fluctu- 
ation transfer function comes from jL,t- Under the con- 
dition V,n <S^V < Vcr, which is easy to satisfy at large A, 
the first integral in Eq. ( 80b ) is much less than the sec- 
ond one, so that the contribution to the noise from the 
region before the virtual cathode may be neglected. Fur- 



thermore, at sufhciently high A, « ^^/K 
Thus, one can write 



1 



-dri 



l/423/2y3/4 



V^- VI) V 



-1/2 



(87) 



where we have taken into account that the main contri- 
bution comes at the upper integration limit and made use 
of the asymptotic expansion of the function hy given by 
Eq. (^) . It is also assumed here that for any fixed energy 
the bias is high, V ^ e. It is justified since the range 
of valuable energies is limited by the Maxwellian expo- 
nentially decaying distribution. Now, substituting the 
Langmuir expression (|5^) for the current and neglecting 
Xm and Vm, one obtains 



7L,t(e) 



'e - 



This formula, after the integration ewer the energies, 
leads to North's asymptotic formulaB for the noise- 
reduction factor: 



F.ifl-^ 
V \ 4 



1.9314 



V 



(89) 



This formula is universal in the sense that it is free 
from any diode parameter including the screening pa- 
rameter A. However, it is assumed that A should be 



sufficiently high to satisfy the simultaneous conditions 
— > oo and V < Vcr- As it is seen from Fig. ^, the 
noise-reduction factor F approaches this asymptotic for- 
mula at high values of the parameters: A > 10^, V > 10"^. 
As we have already noted, in semiconductors it is hard 
to maintain the ballistic regime at biases V ^ 50 because 
of the increasing significance of electron-phonon interac- 
tions which destroy the ballistic regime. In the range of 
interest I ^ V ^ 50 the noise level is seen to be signifi- 
cantly lower than North's asymptotic curve. This means 
that the full set of formulas are necessary to describe 
properly the noise intensity in the semiconductor ballis- 
tic diodes. Another important conclusion from Fig. ^ is 
that for a nondegenerate electron gas there exists the low- 
est noise-reduction level dependent only on the bias and 
the temperature through the factor (/[//(fc^T), and it is 
impossible to surmount it by any choice of the material 
parameter and/or geometrical parameters of the diode. 
This universal minimal-noise curve approaches North's 
asymptotic curve at high voltages. 



H. Spectroscopy of shot noise 

A great advantage of the derived formula (|7^) for the 
current-noise spectral density is that one may obtain the 
partial contribution to the noise from electrons of differ- 
ent injection energies by computing the current fluctu- 
ation transfer functions Jki^)- The electrons for which 
7fe(e) < reduce the current fluctuations. For instance, 
the right-contact electrons always reduce them, since 
7fl(e) < 0, Ve. The reflected carriers originated from 
the left contact (e < Vm) also provide negative values for 
the transfer function and compensate the current fluc- 
tuations by virtue of the potential-barrier fluctuations. 
The same effect is produced by the left-contact transmit- 
ted electrons with the energies slightly above the bar- 
rier height Vm- From both groups, the most efficient 
compensation carriers are those_with the energies in the 
vicinity of Vm where 7l —> — ooH They provide an over- 
compensation of the injected from the contacts fluctua- 
tions. In contrast, the injected electrons whose energy 
greatly exceeds Vm produce negligible perturbations of 
the potential barrier, thus leading to the asymptotic be- 
havior 7L(e) 1, 7_R(e) ^ — 1 as e ^ oo. There also 
exists the specific energy e* , for which the compensation 
fluctuation is exactly equal to the injected fluctuation, 
giving no noise at all, 7l(£*)=0. This curious fact is il- 
lustrated in Fig. |ll| where we present the contribution to 
the current-noise spectral density from different energies 
of electrons injected from the left contact by plotting the 
fimction ^\{e)e~'^ - At high biases, just after the peak 
at £ = Vm, the point with zero contribution to the noise 
is observed. While at equilibrium the maximum contri- 
bution comes from the carriers injected with e = Vm, 
at high biases, when the noise reduction is significant, 
the main contribution comes from the electrons that are 
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injected above the potential barrier height by the value 
about ksT. Therefore, the integral noise-reduction effect 
is a consequence of the suppression of the contributions 
from the electron energies in the vicinity of e*. 
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FIG. 11. Partial contributions to the current-noise spec- 
tral density from different energies e of electrons injected from 
the left contact for biases ^=0.01; 10; 50. The results for 
the right-contact electrons are approximately the same for 
V=0.01 and negligibly small for V^=10 and 50. 



The obtained exact solutions allows us to investigate 
in great detail the correlations between different groups 
of carriers. While the injected carriers are uncorrelated, 
those in the volume of the conductor are strongly cor- 
related, as follows from the derived formulas for the 
fluctuation of t he dist ribut ion function [see general ex- 
pressions ( A12 ) and ( A13 ) in the Appendix]. Those 
correlations may be observed experimentally by mak- 
ing use of a combination of two already realized tech- 
niques: a hot-electron spectrometerEj'Eil and shot-noise 
measurements. EjEj The electron spectrometer, placed 
behind the receiving semitransparent contact, acts-.as.an 
analyzer of electron distribution over the energy.EHlEil In 
this way spectroscopic information, that is, the average 
partial currents /(e) and their fluctuations 61 (e), may be 
measured for different energies e of electrons collected at 
the contact. This is similar to the energy-resolved noise 
measurements realized in Ref. The partial current of 
the transmitted electrons at the receiving (right) contact 
is given by /(e) = /ce^'^^^"6'(e), where the threshold 
energy e=0 corresponds to the arriving electrons that 
have a zero longitudinal kinetic energy at the potential 
minimum. To find the fluctuation (5/(e), we consider 
the fluctuation of the distribution function 6f{x,w) at 
X — X. Since S'ip{X)=0, the terms with 5il)(x) vanish. 
Thus, for the transmitted over the barrier electrons which 
contribute to the current, from Eqs. ( |Al2a| ) and ( |Al3a| ) 
one obtains 



'5/L,i(A, w) = 5fL{X, w) e{w - wr) 

-^e~^"+'' 5V„^:^Siw-WR). (90) 

Since only the positive velocities are considered, one can 
change the velocity variable to the energy by e = uP'—wj^, 
and obtain 

6fLAe) = 5fL{e + V,n) 0{e) - ^e"^^"^" SVmSie). (91) 

By using the relation ( |5^ ) between the fluctuation of the 
contact distribution function and that of the contact in- 
jection current, we obtain 



6Iie) = SlLie + Vra)9{e)-Ic 



5V^5{e). (92) 



Thus, the correlator for the current fluctuations becomes 



(<5/(e)«(e'))U=A-(<5/L(e 



-/, 
-he 



V,n)5lL{e' + Vrn)) 
-''-S{e'){SlLie + Vra)6Vrr^ 
SVrr^ 



-^'"<5(e) {6lL{e' + Vr 
-2^'"<5(e)<5(e') {SV,l) 



(93) 



where the average is taken over the injected fluctuations. 
It is clear that for e, e' > the carriers remain uncor- 
related since only the first term does not vanish. It is 
oc i5(e — e') due to the imposed injection conditions that 
should lead to the full shot noise. In such a case, an 
interesting question arises: What is the reason for the 
noise reduction obtained for the total (integrated over 
the energies) current fluctuations? The answer is found 
looking at the electrons with energies close to the thresh- 
old energy e=0 ("tangent" electrons). All other electrons 
are anticorrelated with that group. This means that if 
there is a positive fluctuation of overbarrier electrons, 
there should be a negative one for the "tangent" elec- 
trons and vice versa. This anticorrelation explains the 
overall noise reduction. The tangent electrons can be 
thought as overcorrelated. The dispersion {6P{i)) has 
a sharp peak at e=0 and then decreases with energy at 
e > 0. This peak is divergent ((5-shaped) in our collision- 
less theory. A small probability of scattering will lead to 
its broadening and finite magnitude. Therefore, by mea- 
suring the dispersion of the partial current fluctuations 
and/or their cross-correlations, one may observe a sharp 
peak and an anticorrelation of electrons, thus making the 
Coulomb correlations effect visible. 



VI. SUMMARY 

In conclusion, we have presented a self-consistent the- 
ory of electron transport and noise in a ballistic two- 
terminal conductor under the conditions of nondegener- 
ate electron gas. Our description is valid for ballistic elec- 
trons in solids as well as in vacuum. By solving analyti- 
cally the kinetic equation coupled self-consistently with a 
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Poisson equation, we have derived the electron distribu- 
tion function and its fluctuation at arbitrary section x of 
the conductor. This allowed us to obtain the steady-state 
spatial distributions of the transport characteristics, the 
I-V curves, and the noise characteristics. While the time- 
averaged quantities are not affected by the Coulomb cor- 
relations, the noise characteristics are demonstrated to 
be drastically modified when those correlations are taken 
into account. Our results are in excellen±.|a£reement with 
the preceding Monte Carlo simulations.Bo 

The obtained formulas have been analyzed in a wide 
range of biases and compared with the correspondent the- 
ory for the vacuum diode. In particular, we have demon- 
strated that the known formulas for vacuum electronics, 
such as the Child 3/2-power law for I-V characteristics 
or North's asymptotic formula for the noise may not be 
applied for the semiconductor diode at biases that are 
relevant for the ballistic transport regime. Instead, one 
should use the more general formulas described in the 
present paper from which follows (i) the linear or sub- 
linear I-V dependences even under a strong limitation 
of transport by a space charge; (ii) a noise level signifi- 
cantly below the level obtained from the North's formula; 
(iii) the sharp recovering of the full shot-noise level at a 
certain critical voltage; (iv) all the transport and noise 
characteristics determined by two dimensionless param- 
eters: the screening parameter A — d/L'jj and the bias 
qU/kBT. 

Based on the derived formula for the current-noise 
spectral density one may distinguish the relative contri- 
butions to the noise from different groups of carriers. We 
have analyzed the contributions coming from the trans- 
mitted and reflected groups of carriers, as well as the 
partial contributions to the noise from electrons of dif- 
ferent injection energies. Finally, it should be noted that 
the analytical approach that we have presented in the 
paper may be extended and applied to various systenia, 
e.g., for different kind of statistics of injecting electrons,ca 
and other types of the contacts. On the other hand, the 
shot-noise suppression effect, which we treat analytically, 
may lead to important applications for low-noise small- 
size semiconductor devices, generators of sub-Poissonian 
light sources,E£l etc. Our work then offers new perspec- 
tives on the study of Coulomb interactions and noise in 
small-size ballistic devices, such as ballistic transistors, 
point contacts, etc. 



APPENDIX A: FLUCTUATIONS OF THE 
DISTRIBUTION FUNCTION AND ELECTRON 
DENSITY 

In a similar way to the subdivision of thestationary dis- 
tribution function (^) into the components correspond- 



ing to different groups of electrons classified in Sec. IV A 
the fluctuation 6f(x,w) may be expressed as 



(Al) 



The boundary conditions for these functions are obtained 
by perturbing the steady-state boundary conditions ( [25| ) 
and using dwk/dipm = — (2iffe)^^. One gets 



SfL.tiO,^c) 

6fL,r{0,Wc) 
5fR.t{\Wc) 



SfbiWc) 0{Wc - Wl) 

+ 7: fL{Wc) S{Wc - wl) Slpn 

5fL{w,)e{wl~wl) 



1 



wr) 

Ir{Wc)5(-Wc - Wr) Siprn, 



2wr 

5fR{w,)e{wl~wl) 

-fR{Wc)S{wR - wl) Slprri- 



(A2a) 

(A2b) 
(A2c) 

(A2d) 



where the additional terms proportional to Sipm describe 
the changes in the distribution functions due to the po- 
tential barrier variation. 

Now we have to solve the perturbed kinetic equation 
(^), which may be rewritten as 



d_ 

dx 



2 dx dw 



5f{x,w) 



1 df dSip 

2 dw dx ' 



(A3) 



where the rhs is supposed to be a given function (for this 
step of calculations). A general solution of this nonhomo- 
geneous partial differential equation is a sum of a solution 
of the homogeneous problem and a particular solution of 
the nonhomogeneous problem. Explicitly, 



CP c rhom I r rnh 



k^L,R, j^t,r. (A4) 



The solution for the homogeneo us p roblem is deter- 
mined by the boundary conditions ([A^). By making use 
of the energy-conservation law (p7|), we make a replace- 
ment 
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Wc = sgn(w) \Jw'^ - Tpix) -I- 



(A5) 



and obtain different contributions to Sfj!"™ in the form 

5/,^°™(x, w) = Sfkix, w)e{±w- w^x)) 

±-^fk{x,w)S(±w - w^{x)) Silj,n, (A6a) 
2w 

6fj:°r{x, w) = Sfkix, w) 9{wUx) - w^) 

-fkix,w)S(wl{x) ~w^)64!m, (A6b) 
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where 5f^°^ and Sf^°^ are defined in the regions < 
X < Xm and Xm < X < X, respectively. The upper sign 
applies for Sf^°^ and the lower sign apphes for Sf'jfJ^, 
both terms vahd in the whole range < a: < A. The 
critical velocity w*(a;) is given by Eq. (p9|). 

The solution of the nonhomogcneous problem can eas- 
ily be found through the steady-state distribution func- 
tion /(et) = f{w^ ~ in terms of the total energy 
Et or, equivalently, in terms of the injection velocity Wc, 



5fZ 



'^'^ dst 
Differentiating Eqs. (]25|), we find 



1 

2w, 



1 S^^--^^54>. 



dWr 



(A7) 



nhojn 
L,t 



5f\ 



nhom 
R,t 



fL{Wc)5'4>{x) 

X 9{wc- wl) - -^5{wc - wl) 
2wc 

fniwc) 6il;{x) 

X 9{-Wc - Wr) + -;:^S{-Wc- Wr) 



2wc 



(A8a) 



(A8b) 



x[9{wl^wl)+6iwl^wl)]. (A8c) 

In these equations the substitution ( [A5| ) is assumed, so 
that the fluctuations are finally the functions of {Xjw). 
Notice that the components for the reflected groups of 
carriers are defined in the regions: 6f2^°"^ for < x < 
Xm and dj"^"™ for a;™ < x < A, while those for the trans- 
mitted groups of carriers are given in the whole range 
< x < A. 

According to the electrostatic boundary conditions 
( p7| ) the fluctuations of the potential at the contacts are 
equal to zero, which leads to vanishing contributions (AS) 
at the contacts (5/"''°"(0, = ^/"''°"(A, w) = 0. The 
contributions (A6) satisfy the boundary conditions ([A^). 
Thus, the distributi on functi on i n the form (A4) with 
eight contributions (AC) and ( |A8| ) is the solution of the 
problem for a given electrostatic potential "^(x) -\- Sip{x). 

For convenience of further consideration, we present Sf 
as a sum of the "injected" and "induced" contributions 



r rind 
"Ik 



L, R, j = t, 



(A9) 



Ik,] "Jk,j : 

In terms of the contact velocities Wc (presented in such a 
form these equations will be frequently used throughout 
of the paper), those contributions are given by 

Sfl7M - Sfkiwc) 0{±wc~ wk), (AlOa) 
Sfk^i^c) = 5h{w,) e{wl - wl), (AlOb) 

and 

5flf{x, w,) = fkiwc) {e{±w, - Wk)54,{x) 

TT^5{±Wc-Wk)[5il}{x)-6tl)j\ , (Alia) 
2wc J 

5f^,\x,w,) = h(w,) {e{wl - wl)5^{x) 

+5{wl - wl)[Si;{x) - • (AUb) 



where the substitution ( [A5| ) is assumed. The same terms 
as functions of (x, w) are determined by the formulas 

Sfiyix,w) =Sfkix,w)e{±w-w4x)), (Al2a) 
Sfiy{x,w)=6fk{x,w)9{wl{x)-w'), (A12b) 

and 

Sfrfix,w) = ^e-'"'+'^(")-''"= {9{±w-w4x))S^Pix) 

Ti^5{ ±w- w4x))[d4>{x) - <5Vm]| , (Al3a) 
2w J 

Sfr:'{x,w) = 1 {9{wl{x) - w')Sij{x) 

+5(wl{x) ~ w^)[6^{x) - Si^m]} ■ (A13b) 

Apparently, has a meaning of the distribution func- 

tion of randomly injected electrons, while Sf^""^ describes 
the change in the steady-state distribution induced by 
injected electrons. 

The obtained fluctuations of the distribution function 
allows one to compute each contribution to the fluctu- 
ations of the electron density Sn{x) by integrating over 
velocities. Changing the integration over w to that over 
the contact injection velocities Wc, we find 



Suk^jix) 



Sfk,jix,w)dw 



Sfk,j {Wc)WcdWc 



oo sgn(wc) \/w'^ + i^ix) - ij}k 



(A14) 



Thus, by using Eqs. ( AlO) and (All), one obtains for the 
injected density fluctuations 



Sn^^Ax) 



oo 



5fL{Wc)Wc dWc 



^^'^ SfR{Wc)WcdWc 



Sfk{Wc)WcdWc 



(Al5a) 
(A15b) 
(A15c) 



and for the induced fluctuations 

^nT^ix) = nk,t S^pix) - — — 



■[Slp{x) - Sipr, 



Sn'j^^ix) = nk,rS'ip{x) 



f^TW.^,{x) 



[5il){x) - 6ip„ 



(Al6a) 
(A16b) 



Here, the contributions (A15) can be interpreted as the 
electron-density fluctuations at a slice x caused by the 
stochastic injection from the contacts to the base. The 
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contributions ( A16) are related to a variation of the sta- 
tionary electron density due to a local variation of the po- 
tential and its minimal value (a self-consistent response) . 
As before, the terms SriL^r and Snu^r are defined on the 
intervals < a; < Xm and Xm < x < X, respectively, while 
the terms Snk^t are defined on the whole range < a; < A. 
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